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ABSTRACT 


A  model  that  produces  three-dimensional  ray  trajectories  in  an  optical  fiber  is 
derived  through  the  use  of  a  simple  correspondence  rule.  The  analysis  and  simulations 
presented  will  be  in  terms  of  dimensionless  parameters.  The  curves  prepared  are 
presented  for  allowed  radial  and  azimuthal  mode  numbers  and  are  presented  in  terms  of 
dimensions  which  are  normalized  by  the  core  radius.  The  model  presented  will  be  shown 
to  be  in  agreement  with  the  standard  Eikonal  ray  theory  and  will  serve  as  a  link  between 
the  ray  trajectories  and  the  mode  numbers  of  the  WKB  waveguide  solutions  of  an  optical 
fiber. 
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I. 


INTRODUCTION 


A.  REVIEW  OF  PRIOR  WORK 

An  extensive  body  of  literature  exists  for  both  wave  and 
ray  characterization  of  optical  fibers.  It  is  generally 
accepted  that  exact  waveguide  solutions  exist  for  a  select 
number  of  specific  index  profiles  [Ref.  l,  2]  .  The  exact 
waveguide  solutions  are  frequently  derived  for  step  index 
fibers  in  intermediate  level  optical  fiber  texts.  For  example, 
Cherin's  Optical  Fiber  Communications  [Ref.  3]  provides  a 
fairly  clear  waveguide  analysis  of  the  step  index  fiber.  In 
the  general  case  of  graded  index  fibers  the  problem  is  not 
analytically  tractable  without  approximation  [Ref.  1].  The 
most  common  technique  for  the  analysis  of  graded  index  fibers 
is  based  on  the  WKB  approximation  [Ref.  4,  5,  6,  7]  .  It  is 
possible  to  show  that  the  first-order  WKB  approximation  is 
equivalent  to  the  results  produced  from  a  simple  ray  model 
[Ref.  8]  . 

Despite  the  mathematical  importance  of  having  an  exact 
solution  for  the  waveguide  fibers  in  an  optical  fiber,  the 
propagation  characteristics  can  be  difficult  to  visualize.  It 
is  not  surprising  that  an  alternative  approach  to  treating 
optical  fibers,  based  on  the  Eikonal  ray  approach  [Ref.  9,  10, 
11,  12]  ,  which  in  turn  is  based  on  Fermat's  extremum  principle 
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[Ref.  10,  11,  12,  13,  14],  has  proven  to  be  extremely  useful 
in  explaining  effects  which  take  place  in  optical  fibers.  For 
example,  application  of  such  techniques  to  the  prediction  of 
bending  losses  in  optical  fibers  has  appeared  in  the 
literature  [Ref.  13,  15]  .  In  addition,  achromatic  modal 
dispersion  effects  have  been  predicted  with  ray  models  [Ref. 
16]  .  It  has  been  shown  that  the  exact  step  index  radical 
function  has  a  first  inflection  point  at  the  inner  ray  caustic 
or  turning  point  [Ref.  8].  To  a  significant  degree  the 
association  between  the  WKB  wave  solutions,  which  depends  on 
azimuthal  and  radial  mode  numbers,  and  ray  models  has  not  been 
fully  explored.  For  example,  the  ray  trajectories  were 
typically  derived  for  arbitrary  initial  conditions  without 
regard  for  ray  conditions  dictated  by  invariant  conditions 
associated  with  mode  specific  quantum  numbers. 

Recently,  an  approach  based  on  a  simple  correspondence 
between  waves  and  rays  has  been  proposed  [Ref.  17]  as  a 
substitute  for  the  first-order  WKB  approach  of  counting  modes 
[Ref.  9]  .  This  correspondence  principle  has  been  used  to 
qualitatively  demonstrate  that  there  should  be  a  direct 
association  between  the  mode  numbers  in  an  optical  fiber  and 
the  ray  trajectories.  The  main  focus  of  this  thesis  is  to 
treat  this  association  quantitatively. 
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B.  OBJECTIVES  OF  THE  THESIS 

The  correspondence  principle  [Ref.  17]  leads  to  rules 
which  predict  both  initial  conditions  and  propagator 
conditions  on  ray  trajectories  associated  with  specific  modes. 
This  thesis  will  cover  all  background  analysis  leading  to 
application  of  these  rules  to  meridional  rays  and  skew  rays  in 
a  graded  index  fiber.  Whenever  necessary  for  visualization, 
full  three-dimensional  ray  trajectories  are  provided.  The 
analysis  and  program  simulations  are  derived  and  implemented, 
respectively,  in  a  normalized  form  with  respect  to  the  vacuum 
wavelength  and  core  radius.  As  a  result,  the  computer  results 
are  generalized  to  a  broader  class  of  problems.  Specifically, 
the  optical  fiber  rays  are  expressed  in  terms  of  the  commonly 
used  V-paraimeter,  numerical  aperture  and  core  index.  Appendix 
B  has  been  dedicated  to  demonstrating  that  the  correspondence 
principle,  mentioned  cibove,  leads  to  propagator  conditions 
which  are  precisely  equivalent  to  rules  obtained  from  the 
standard  Eikonal  model.  Appendix  A  has  outlined  in  flowchart 
fashion  the  conceptional  development  of  the  programs.  Appendix 
E  of  the  thesis  manuscript  has  been  scheduled  for  an  IEEE 
publication  [Ref.  19]  in  March  1993. 

C.  OVERVIEW  OF  THE  THESIS 

Chapter  I  is  the  introduction  and  covers  both  a  review  of 
relevant  literature  and  a  thesis  overview.  In  the  interest  of 
keeping  this  thesis  fairly  self-contained.  Chapter  II  provides 
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a  brief  technical  description  of  essential  background 
preliminaries . 

In  Chapter  III,  the  basis  for  the  specific  problems 
examined  in  later  chapters  is  provided.  A  conceptual  flowchart 
shown  in  Appendix  A  summarizes  the  basic  analysis  which  leads 
to  mode  specific  trajectories. 

In  Chapter  IV,  attention  is  given  to  meridional  rays  for 
the  arbitrary  a  profile.  The  specific  meridional  ray  cases  of 
step  index  fibers  (q!=<»)  and  parabolicity  (q!=2)  are  both 
analyzed  in  detail.  Appendix  D  provides  a  useful  mathematical 
analysis  for  the  essential  function  F(a),  a  normalization 
integral  used  in  this  chapter. 

In  Chapter  V,  the  solutions  for  the  step  index  skew  rays 
are  obtained.  This  chapter  is  a  link  between  the  general 
formulation  of  Chapter  III  and  the  more  complete  step  index 
analysis  contained  in  Appendix  E. 

In  Chapter  VI,  numerical  solutions  for  skew  rays  with 
arbitrary  of  profile  is  provided.  Some  useful  variations  of 
fiber  optic  parameters  due  to  initial  conditions  are  also 
demonstrated  and  explained  here. 

In  Chapter  VII,  quasi  analytic  solutions  for  the 
interesting  case  of  skew  rays  with  a=2  are  presented.  The  well 
known  fact  of  the  orbital  periodicity  of  this  case  is  also 
proven  here. 

Appendix  B  shows  the  exact  agreement  between  the  ray 
analysis  presented  and  the  Eikonal  analysis. 
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Appendix  C  provides  a  symbol  table  which  contains  notation 
for  the  most  common  symbols  used  throughout  in  the  thesis. 
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II 


BACKGROUND 


Starting  from  the  Hemholtz  equation,  which  is  the  phasor 
form  of  the  linear  wave  equation  for  isotropic,  source-free, 
homogeneous  media,  we  have 

[V2+Jc2(r)  ]B=0,  (2-1) 

where  k(r)  is  the  wave  number  and  E  is  the  electric  field.  It 
is  noted  that  this  equation  is  only  approximately  valid  in  a 
graded  index  fiber.  Consistent  with  the  form  given  in  (2-1) , 
we  assiome  that 

k  ( r )  =n  ( r )  +k^e^,  (2-2) 

where  k  is  the  vacuum  wave  ntutiber,  n{f)  is  the  index  of 
refraction,  (r,  0,  z)  are  the  cylindrical  coordinates  needed 
to  describe  points  in  the  fiber  and  (e,.,  e^)  are  the  unit 

vectors  of  the  coordinate  system. 

We  list  here  the  basic  assvimptions  that  n(f)  satisfies 

1.  The  index  depends  on  the  distance  from  the  symmetry  axis 
of  the  fiber  (viz.  n(r)=n(r)). 

2.  The  index  is  a  decreasing  function  for  rsa. 

3 .  The  index  is  constant  for  raa . 

These  restrictions  can  be  represented  in  Figure  l  where 
k(r)=kn(r)  . 
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Figure  1:  Restrictions  for  k(r) . 

Following  a  standard  procedure  [Ref.  l]  in  waveguide 
analysis,  we  assume  that  the  general  phasor  for  the  z 
component  of  the  electric  field  in  the  cylindrical  optical 
fiber  can  be  represented  as 

E,(r)  =F,(r)  (2-3) 

where  F,(r)  is  the  radial  solution  to  (2-1),  /S  is  the  phase 
constant  of  the  waveguide,  and  v  is  the  azimuthal  mode  number, 
which  is  forced  to  be  an  integer  due  to  periodic  boundary 
conditions  [Ref.  l] . 
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In  the  case  of  a  step  index  fiber,  an  application  of 
Maxwell's  equations  leads  to  a  transcendental  equation  for  the 
phase  constant.  The  equation  has  a  finite  number  of  discrete 
solutions  for  each  integer  u  [Ref.  l]  .  The  obtained  convention 
for  ordering  these  solutions  is 

/3  >/3, . . .  >0,  >i8,  ,  (2-4) 

I  a  m 

where  the  second  subscript  m  is  known  as  the  radial  mode 
number.  This  convention  is  also  applied  to  the  approximate  WKB 
[Ref.  7]  solutions  in  inhomogeneous  optical  fibers. 

A  simple  correspondence  rule  for  plane  waves  [Ref.  17] 

v--jk,  (2-5) 

can  be  applied  to  the  <f>  and  z  dependencies  of  (2-3)  .  The  grad 
operator  is  given  by  [Ref.  18] 


"iff— a  ^  ^  j.a  ^ 


It  now  follows  from  (2-6)  and  (2-3)  that 


(2-6) 


J 


which  from  (2-5)  leads  to 


(2-7) 


Similarly, 


j{vE,)^=-j0E^. 


and  again  from  (2-5) 


(2-8) 


(2-9) 
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(2-10) 

It  is  noted  from  (2-3)  that  the  correspondence  rule  applies 
only  to  the  4>  and  z  components. 

Using  the  relation  between  the  propagation  vector  and  its 
components,  it  follows  from  (2-2),  (2-8)  and  (2-10)  that  the 

radial  component  of  the  propagation  vector  can  be  expressed  as 


kHr)  . 

r 


(2-11) 


Both  +  and  -  roots  are  required  for  the  development  of  a  full 
cycle  of  a  trajectory.  Nevertheless,  symmetry  conditions  in 
the  trajectory  will  allow  the  analysis  to  be  based  on  the  + 
root  only.  Furthermore,  for  k,  to  correspond  to  a  bounded 
orbit,  it  must  be  real  [Ref.  l,  8] .  Without  loss  in 
generality,  we  can  continue  with  only  the  +  form  for  k,. 
According  to  (2-4)  and  (2-11),  the  radial  component  of  the 
propagation  vector  will  increase  with  increasing  m  (see  Figure 
2) .  The  crosshatched  domain  of  Figure  3  represents  the 
condition  that  k,  is  real.  Therefore,  terms  under  the  radical 
of  (2-11)  satisfy 

kr=k^(r) -0^- {  —  )‘^z0 .  (2-12) 

A  condition  of  phase  synchronization  [Ref.  1]  on  the 
radial  part  of  the  wave  vector  results  in  the  approximate 
expression 
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Figure  2;  Radial  component  of  propagation  vector  k,  versus  m. 


k^(r) dr,  (2-13) 

r 

where  the  radial  mode  number  m  must  take  on  positive  integer 
values.  The  integral  limits  r,  and  r2,  are  the  two  values  for 
which  the  integrand  of  (2-13)  is  zero.  These  are  generally 
referred  to  as  the  turning  points  as  represented  in  Figure  3 . 
We  note  that  r,  and  rj  satisfy 

k^=kHr)  -0^- (1)^=0.  (2-14) 

r 

More  precisely,  the  relation  (2-13)  becomes  exact  after  the 

substitutions  m*ni+l  and  [Ref.  5]  .  In  the  WKB  analysis 

2  4 

of  optical  fibers,  these  offsets  are  usually  ignored. 
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A  significant  point  relates  Figure  3  to  the  radial  mode 
nvimber  given  by  (2-13)  .  The  value  of  m  is  qualitatively 
proportional  to  the  cross-hatched  area.  The  general  procedure 
to  find  the  maximum  m  is  to  allow  to  be  equal  to  If  /3^ 
falls  below  then  the  additional  modes  created  are  no 
longer  trapped.  Many  authors  have  explored  the  interpretation 
that  the  optical  fiber  trajectory  is  equivalent  to  the 
trajectory  of  a  particle  in  a  potential  well  [Ref.  8] . 

Figure  3  correctly  suggests  that  there  is  a  maximum  p  for 
which  the  cross-hatched  region  vanishes.  Clearly  there  is  an 
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upper  limit  on  the  value  of  v  which  occurs  for  the  condition 
r,=r2.  The  general  relation  between  the  turning  points  and  p 
can  be  obtained  from  (2-14)  as 


v=z\lk}  (r) ,  (2-15) 

which  is  represented  in  Figure  4.  It  is  apparent  from  Figure 
3  that  the  condition  fi=kn2  will  allow  the  maximum  range  in  p 
and  the  corresponding  maximum  in  the  peak  on  Figure  4 .  In 
calculating  p^,  p  (r^)  will  in  general  be  evaluated  under  the 
stated  limiting  condition  /3=kn2.  At  p=Q ,  the  solutions  for 
turning  points  are  easily  read  off  in  Figure  4  as  0  and 

n'M-^).  Since  J?=n,cos0,  the  second  condition  leads  to  the 
k  k 

condition  for  total  internal  reflection 

n(r2)  =n,cosP.  (2-16) 

For  the  step  index  fiber,  r2=a,  which  leads  to  the  well  Icnown 
relation 


003$  =  —,  (2-17) 
n, 

Equations  (2-8),  (2-10),  (2-11),  (2-12)  and  (2-15)  become 
the  basis  for  the  three-dimensional  ray  modeling  that  will  be 
the  focus  of  the  discussion  that  follows.  It  is  noted  that  the 
components  (k,,  k^,  )q)  will  describe  the  direction  of  the  ray 
and  that  coordinates  (r,  0,  z)  locate  the  ray.  For  purposes  of 


Fi^9er4: 4Gle£tez]a}:a±el:aft44}iLahebs£snaetiutniii9ngo{2atLBtandnd.  v . 


anaii$[iys  thBheregBeinetcdiralrelaftdBiiBh^iiapararsioTOreacafsIii^ly 
de  flMededn  itih^haex^xthsp&pter . 


Best  Available 


III.  GENERAL  CHARACTERIZATION  OP  THE  THREE-DIMENSIONAL  RAY 
MODEL 

A.  INTRODUCTORY  REMARKS 

This  chapter  will  provide  a  basis  for  the  specific 
categories  of  problems  examined  in  later  chapters  (e.g., 
meridional,  skew  step  index,  and  skew  graded  index)  .  A  general 
approach  to  the  goal  of  obtaining  mode -connected  ray 
trajectories  will  be  presented  at  the  end  of  the  chapter. 


B.  GEOMETRY  OF  THE  PROBLEM 

With  reference  to  Figure  5,  the  orientation  of  the  wave 
vector  is  defined  in  terms  of  6  and  It  follows  that 


k^=k{r)  8±n$  sin$, 

(3-1) 

kj=k(r)  coad , 

(3-2) 

and 

k^=k{r)  aind  cos{. 

(3-3) 

The  incremental  line  variation  in  the  ray 

defined  as 

trajectory  can  be 

dl=dz  Cj+dr  e^+r  d<^  e^. 

(3-4) 

The  propagation  vector,  defined  by  the  cylindrical  components 

(3-1),  (3-2)  and  (3-3),  is  perpendicular  to  the  wavefronts  of 

constant  phase  and  therefore  tangent  to  the  ray  trajectories. 
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Figure  5:  Geometry  of  the  problem. 


15 


Therefore,  the  incremental  components  d</),  dr,  and  dz  are 
tangential  to  (2-8),  (2-10)  and  (2-11),  respectively,  and  it 
follows 

dz-dl  cos6 ,  (3-5) 

dr=dl  sinO  cos|  (3-6) 

and 

r  d</)=dl  3±n6  sin^.  (3-7) 

C.  MODELING  ANALYSIS 

1.  Computation  of  the  tuxming  points 

Without  any  loss  in  generality,  we  assume  that  the  index 
profile  can  be  expressed  in  terms  of  a  standard  power  profile, 
which  is  matched  to  the  edge  conditions, 

n  ( 0 )  =n,  ( 3  -  8 ) 

and 

n(a)  =nj,  (3-9) 

and  a  perturbation  f(r)  as  shown  below 

kWir)  =k^[n^{r) +NA^f{r)]  ,  (3-10) 

where 

f(0)=f(a)=0.  (3-11) 

The  numerical  aperture  of  the  fiber  is  given  by  [Ref.  1] 
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NA 


and  the  power  profile  satisfies 


(3-12) 


np'=nf[l-2A(£)“] 
^  a 


(3-13) 


in  which  the  index  difference  for  this  formula  is  defined  as 


A= 


2  2 
^1  --02 

2nf 


(3-14) 


It  can  be  easily  checked  that  at  r=a,  np=n2. 

Substitution  of  (3-12)  into  (2-12)  leads  to 


m= 


€(/3)P^f(r)-(-^)^-lF“*^^ 

TrJ;.J  V  r 


e(i8)  = 


NA^ 


and 


(3-15) 


where  the  V-parameter  is  given  by  [Ref.  1] 

V=kNAa,  (3-16) 

the  normalized  phase  constant,  e  (j8)  ,  is  defined  by 


(3-17) 


;3-18) 


defines  the  radial  normalization. 
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In  order  to  further  analyze  the  problem,  the  perturbation 
of  the  index  power  will  be  dropped.  After  taking  f(r)=0,  the 
condition  on  the  turning  points  will  satisfy 

r“*2-€  {/3)P+(Z)2=0,  (3-19) 

which  is  a  normalized  version  of  condition  (2-14)  .  Figure  6  is 
plotted  for  ot=2,  V=100,  *'=10  and  e  (jS)  as  stated  on  the  Figure. 
It  is  noted  that  as  e  (/3)  increases,  the  zero  crossing  points 
(or  turning  points)  should  move  closer  together.  Eventually 
the  turning  points  will  coalesce. 

2 .  Computation  of  the  maximum  radial  mode  number 

As  described  in  Chapter  II,  the  maximum  radial  mode  number 
is  obtained  by  setting  or  according  to  (3-17), 

Thus,  (3-15)  leads  to 

f . 

where  the  turning  points  which  correspond  to  the  maximum  m 
will  then  satisfy 

r '“*^-r  ( ^)  ^=0 .  (3-21) 

3 .  Confutation  of  the  maximum  azimuthal  mode  number 

The  first  step  in  obtaining  the  maximum  azimuthal  mode 
number  is  to  find  i'(rp)  as  represented  in  Figure  4.  There 
is  a  peak  in  the  function  i»(r)  already  defined  by  (2-15)  .  The 
second  step  is  to  evaluate  the  expression  for  the  limiting 
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Figure  6:  Left  hand  side  of  Equation  (3-19). 
condition  /J^knj. 

Combining  the  power  profile  (3-13)  with  (2-15)  leads  to 

[  1 -2A  ( £ )  “]  .  ( 3  -  22 ) 

a 

By  evaluating 

^{p^)=0,  (3-23) 

it  can  be  shown  that 


n. 


n,A(2-t-a) 


(arbitrary  v)  , 


(3-24) 
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which  for  the  limit  condition  ^=kn2,  leads  to 


2  - 
r  =  {  )  “ 


(3-25) 


where  r  is  the  normalized  radius  extremum  given  by 


(3-26) 


Finally,  from  substitution  of  (3-25)  into  (3-24),  we  have 


V  =V{ 

““  2+0! 


1 

)  “ 


a 


2+a 


(3-27) 


which  shows  that  depends  on  the  a  profile  and  the  V- 

parameter . 

4.  Coi^putation  of  the  evolution  equations 

Through  combinations  of  (3-2),  (3-3),  (3-5)  and  (3-6)  we 

have  that 


dz ^dz 0  _  ctnd 

dr  dr  K  K  cos^ ' 


where 


z=£  (3-29) 

a 

is  the  normalized  z  component  of  the  ray. 

It  is  easy  to  demonstrate  from  (2-11),  (2-12)  and  (3-15) 

that 
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(3-30) 


e  W)  r^+f{r)  -  ( p  , 

and  from  (3-16)  and  (3-17)  that 

ff_  \lni-NA^e  ip)  (3-31) 

^  a  NA 

Substitution  of  (3-30)  and  (3-31)  into  (3-28)  leads  to 

(iz_  \ln^-NA^€  (/S)  r 

Wft  I - - - ~ '  (3-32) 

e  (|S)r2-(^)2-(r)2- 

which  is  the  evolution  equation  for  the  normalized  z  component 
of  the  ray. 

Wor)cing  similarly  for  the  normalized  <{>  component  of  the 
ray,  combination  of  (3-1),  (3-3),  (3-6)  and  (3-7)  gives 


dr 


V 

r 


=tan^ , 


and  after  substitution  of  (3-30)  leads  to 


d4>  _  _  V  1 
dr  V  r 


_ 1 _ 

€{/3)P-(-^)^-?“*^ 

>1 


(3-33) 


(3-34) 


which  is  the  evolution  equation  for  the  normalized  4>  component 
of  the  ray.  Equations  (3-32)  and  (3-34),  which  describe  the 
ray  trajectories,  are  shown  to  be  in  agreement  with  the 
Eikonal  solution  [Ref.  9]  in  Appendix  B.  Unlilce  the  Ei)conal 
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evolution  equations,  (3-32)  and  (3-34)  are  expressed  in  terms 
of  the  azimuthal  mode  number  and  the  phase  constant  associated 
with  the  radial  mode  number  (2-13) . 

The  procedure  describing  the  general  approach  applied  in 
subsequent  chapters  is  summarized  in  Appendix  A. 
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IV.  MERIDIONAL  RAYS 


In  this  chapter  attention  is  given  to  meridional  rays.  It 
is  well  known  that  modal  dispersion  can  be  approximated  fairly 
well  through  consideration  of  the  meridional  rays  only  [Ref. 
1]  .  In  the  first  section,  the  treatment  for  an  arbitrary  a 
parameter  is  presentad.  In  the  last  section,  some  conclusions 
for  the  specific  case  of  step  index  fibers  are  then  presented. 


A.  GRADED  INDEX  MERIDIONAL  RAYS  WITH  ARBITRARY  a 

Meridional  rays  are  rays  which  pass  through  the  z  axis, 
and  according  to  (2-8)  ,  have  no  rotational  component,  that  is, 

i»=(K»k^=0.  (4-1) 

For  a  power  profile  given  by  (3-13)  ,  it  follows  from  (4-1) 
that  (3-15)  and  (3-20)  simplify  for  this  case  to 

jn=  if fVe  ( )8 ) ^  (4-2) 

ttJo  r 

and 


IT  jo 


r 


(4-3) 


respectively.  It  is  noted  that  the  lower  limit  of  integration 
in  (4-2)  and  (4-3)  is  zero  since  the  associated  integrands 
have  a  root  at  r=0 .  This  agrees  with  the  physical 
interpretation  consistent  with  Figure  7,  which  represents  the 
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Figure  7;  Physical  interpretation  for  meridional  rays. 


special  case  i'=0  of  Figure  3.  Equation  (4-2)  can  be  rewritten 
as 


IT 


r 

Jo  J 


1- 


TWf 


dr, 


(4-4) 


and  therefore,  the  second  turning  point  will,  from  (4-4) , 
satisfy 


2 

?2=(€()3)) 

After  the  variable  substitution 


(4-5) 
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u=^  =  . 


T' 


(4-6) 


{€  (jS)  ) 


we  find 


24a 


m 


ir 


(4-7) 


where 


F(of) du. 


(4-8) 


Although  a  general  solution  to  (4-8)  for  all  a  does  not  appear 
possible,  the  integral  F(oi)  is  easily  computed  via  numerical 
simulation  (see  Figure  8).  The  functions  F,  (a)  and  F2(q!) 
defined  in  Tcible  1  are  seen  in  Figure  8  to  agree  with  P(o;)  for 
small  and  large  a,  respectively.  The  analytic  properties  of 
F{(x)  are  discussed  in  detail  in  Appendix  D.  Table  1 
summarizing  this  discussion  is  provided  here. 


a 

F(a) 

1 

2/3 

2 

7r/4 

0 

0  (non-physical) 

00 

1 

QOl 

F,  (cy)  =  (2a+l)/(2a+2) 

0f<l 

F2(a)  =  (a7r)''V2 

Table  1:  Values  for  F(a). 
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lO 


Figure  8:  F(oif) 


In  order  to  evaluate  from  (4-3),  e  (/3)  is  again  set  to 


one  (see  Chapter  III  Section  B) ,  in  (4-7),  leading  to 

(4-9) 

A  representative  meridional  graded  index  ray  trajectory  is 
shown  in  Figure  9. 

Before  deriving  evolution  equations  for  ray  trajectories, 
we  observe  from  (2-10)  and  (3-2)  that 

j8=Jc  n(r)  cos0„(r)  =Jc  n,cos6,„  (4-10) 
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where  the  subscript  "1"  is  to  indicate  evaluation  at  n(r)=n, 
(see  (3-8)),  and  subscript  m  is  to  indicate  a  discrete  angle 
associated  with  a  specific  radial  mode.  A  geometrical 
interpretation  which  relates  the  angle  to  m  can  be  made. 
After  substituting  (4-10)  into  (3-17)  we  get 


€(^)  = 

which  when  substituted  into 
between  and  m  as 


NA^ 

(4-7) 


,  (‘i-lD 

provides  a  relationship 


s±nd 


NA ,  -nm  1  ^ 

'“F  fToT 


(4-12) 


which  is  in  agreement  with  Figure  2. 

For  calculation  of  z(r) ,  we  see  that  (3-32)  simplifies 
under  condition  (4-1)  to 


dz_\/ii?+WA^€  (/3) 
NA 

After  combinations  of  (3-17), 
on  the  integral  form 


_ r _ 

v/e  (/S)?-(?)“'*2> 

(4-10),  and  (4-11), 


(4-13) 


z(r)  takes 


z(r)  = 


n|Cosg,„  rr 
n,sind,„Jo 


dr 


(4-14) 


After  use  of  (4-6) ,  it  is  easily  seen  that 
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z{r)  =r.ctnd,J''-  ^4— .  (4-15) 

v'l-u” 

where  re[0,r2]  .  The  parameter  rj  is  defined  in  terms  of  m 
through  (4-5) ,  (4-11) ,  and  (4-12)  which  leads  to 

r={^^_i_)^  (4-16) 

*  '  V  F{ot) 

Again,  as  in  the  case  of  F(cx),  exact  solutions  for  (4-14)  and 
(4-15)  can  be  found  for  q;=1,  ck=2,  and  a=oo. 

Because  of  condition  (4-1)  ,  the  azimuthal  variation  is 
zero.  Therefore,  there  is  no  need  to  consider  the  4>  variation 
described  by  (3-34)  .  Appendix  A  summarizes  the  process  of 
trajectory  analysis  for  meridional  rays  also. 

Figure  10  illustrates  meridional  rays  plotted  using  (4-14) 
under  conditions  a=6,  NA=0.5,  V=100  and  n,=0.5.  The  curves  are 
distinguished  only  by  the  radial  mode  numbers  m*20  and  m=2.  As 
expected  from  Figure  2,  m  increases  with  k,,  and  from  (2-4), 
0  decreases  with  m.  Examination  of  Figure  6  confirms  that  the 
peak  values  of  the  trajectories  increase  with  m.  The  values 
predicted  from  (4-16),  r2=0.91l4  for  m=20  and  r2=0.5125  for 
m=2,  agree  with  the  curves  shown. 

B.  GRADED  INDEX  MERIDIONAL  RAYS  WITH  a-2 

Equations  (4-9),  (4-12),  (4-15)  and  (4-16)  can  be 

evaluated  for  o(=2.  From  Table  1,  (4-9),  (4-12),  and  (4-16) 

become 
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in  o  »n 
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I 
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Figure  10:  Meridional  ray  trajectories  for  a=6. 
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normolized 


(4-17) 


sind,„=i®[  ii?  (4-18) 


and 


r-\  4ni  (4-19) 

-V  ' 

respectively.  And  lastly,  through  use  of  the  integral  [Ref. 
18] 


*  du 

'  v/W 


;=arcsinu. 


(4-20) 


the  trajectory  equation  (4-15)  can  be  expressed  as 


r=r,sin( 


T^ctne,„ 


)  . 


(4-21) 


Therefore,  for  a'=2,  the  trajectory  shape  is  exactly  sinusoidal 
with  a  normalized  z  period 

A^=27rr2Ctn^,^,  (4-22) 

defined  in  terms  of  the  specific  initial  condition  of  the 
angle  of  the  ray  as  it  crosses  the  normalized  z  axis.  From  (4- 
18)  and  (4-19)  ,  the  dependence  of  the  normalized  period  on 
V,NA  and  m  is  given  by 
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(4-23) 


where 


Y= 


4m 

-V’ 


(4-24) 


and 


X={^)\  (4-25) 

NA 

Figure  11  is  a  plot  of  A„  for  typical  values  of  X. 

Figure  12  shows  two  meridional  rays  plotted  using  (4-22) 
under  conditions  a=2,  NA=0.5,  V=100  and  n,=1.5.  The  curves  are 
only  distinguished  by  the  radial  mode  numbers  m=20  and  m=2.  As 
expected  from  (4-19),  the  peak  value  of  the  trajectories  are 
1-2=0.894  for  m=20  and  r2=0.282  for  m=2 .  Also,  it  can  be  seen 
from  the  X=9  curve  of  Figure  11,  that  the  period  for  m=2  is 
approximately  19  and  the  period  for  m=20  is  18. 


C.  STEP  INDEX  MERIDIONAL  RATS 

By  definition,  for  step  index  meridional  rays,  in  the 
limit  as 

a-*oo,  (4-26) 

which  when  substituted  into  (3-27)  gives 

p^=Int[V]  .  (4-27) 

The  peak  turning  point  can  be  found  by  substituting  (4-1) 
and  (4-26)  into  (3-25)  yielding 
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Figure  11 :  A„  versus  m  for  meridional  rays  with  a*2 . 
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ozimuthol  mode  number 


o 

00 


Figure  12:  Meridional  ray  trajectories  for  a=2 . 
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normolized  z 


(4-28) 


In  order  to  find  the  maximum  radial  mode  number  (4-1)  ,  (4- 
26)  and  (4-28)  are  substituted  into  (3-20)  which  gives 

ir\^=lnt[^].  (4-29) 

The  evolution  equation  for  the  z  component  can  be  derived 
by  substituting  (4-26)  into  (3-28)  yielding 

—  =ctn0.  (4-30) 

dr 

The  still  unknown  angle  6  of  (4-30)  can  be  found  from  (2- 
12) ,  which  for  this  case  takes  the  form 

in=— r*’kn,sin^dr.  (4-31) 

ir  Jo  ‘ 

By  integrating  (4-31)  and  eliminating  a  and  k  through 
combinations  of  (3-12)  and  (3-16) ,  we  have 

g=sin~'  ( ■”*^^)  ,  (4-32) 

Vn, 

Finally,  by  integration  of  (4-30),  we  get 

i(r)=rctn(9,  (4-33) 

where  according  to  (4-28),  re [0,1].  Equations  (4-32)  and  (4- 
33)  are  the  basis  for  the  simulations  described  in  the  next 
paragraph. 

Figure  13  illustrates  step  index  meridional  rays  plotted 
using  (4-33)  under  conditions  NA=0.5,  V=100  and  n,=1.5.  The 
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Figure  13:  Step  index  meridional  ray  trajectories. 
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normolized  z 


curves  are  only  distinguished  by  the  radial  mode  values  m=l 
and  m=20.  Although  it  is  not  apparent  from  Figure  13,  as 
e3q)ected  from  (4-28),  in  both  cases,  rp=l. 

D.  FINAL  REMARKS 

In  this  chapter,  the  general  equations  of  Chapter  III  are 
applied  under  conditions  valid  for  meridional  rays.  An 
analytic  form  for  the  ray  trajectories  was  established.  The 
existence  of  a  direct  connection  between  the  radial  mode 
number  m  and  the  orientation  angle  6  has  been  demonstrated. 
This  is  a  first  sign  that  a  link  between  Eikonal  ray 
trajectories  and  mode  numbers  is  possible.  It  is  noted  that 
although  a  three-dimensional  perspective  of  meridional  ray 
trajectories  could  have  been  demonstrated,  the  fact  that  the 
rotational  component  is  zero  (see  (4-1))  made  it  unnecessary. 
In  the  chapters  that  follow,  the  analysis  is  applied  to  skew 
rays . 
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V. 


SKEW  RJLYS  IN  A  STEP  INDEX  OPTICAL  FIBER 


A.  INTRODUCTORY  REMARKS 

This  Chapter  presents  a  brief  discussion  for  the  analysis 
of  skew  rays  in  a  step  index  fiber.  It  follows  from  (3-13) 
that  the  general  results  from  Chapter  III  can  be  applied  after 
taking  the  limit  of-w.  A  more  detailed  description  of  this 
subject  can  be  found  in  Appendix  E.  This  appendix  is  a 
reproduction  of  a  published  work  [Ref.  19]  based  on  results 
obtained  in  the  early  stage  of  the  thesis.  One  purpose  of  this 
chapter  is  to  provide  a  clear  link  between  the  general 
formulation  of  Chapter  III  and  the  analytical  approach  for 
step  index  optical  fibers  described  in  Appendix  E.  Some 
additional  results  will  also  be  presented.  In  this  chapter  the 
notation  on  the  turning  points  will  be  modified  to  be  in 
agreement  with  Appendix  E.  Specifically,  the  more  suggestive 
labels  of  and  r,^  will  be  used  for  r,  and  rj,  respectively. 

B.  OVERALL  ANALYSIS 

In  this  section  the  solutions  in  the  step  index  case  (in 
the  limit  as  Of-*oo)  are  obtained.  It  follows  from  Figure  14  that 
r„„=l.  This  can  be  confirmed  from  (3-25) 
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r  ^  )  “-►1 . 

QfHw'  2+a 


(5-1) 


It  then  follows  from  (3-15),  after  ta]cing  f(r)=0,  that 


m= 


-3lf 


€(/S)r^-(^)'  dr, 


(5-2) 


where  the  lower  limit  can  be  obtained  from  (3-19) 


•rmm= 

fUBl 


vVe  (|8) 


(5-3) 


After  substitution  of  (5-3)  into  (5-2) ,  we  find  that 
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1=-*) 

r  r 

which  can  be  directly  integrated  using  [Ref.  18] 

r  S—  ■  du^V' u^-a^  -asec’  ( — )  ,  (5-5) 

J  u  a 

which  leads  to 

jn=^[r  ( -sec-‘(^)]  .  (5-6) 

^  J  ^min  ^min 

which  apparently  is  not  appliccible  for  meridional  rays  {u=0) 
since  r^  given  by  (5-3)  for  meridional  rays  is  zero.  This 
case  has  been  covered  in  Chapter  IV,  Section  C. 

As  in  Chapter  III,  Section  B,  the  con5)utation  of 
requires  setting  €(/3)=1.0.  It  follows  from  (5-2)  that  for 

(5-7) 

and  the  corresponding  value  of  m^  is 

iit...  ( y )  =-^  -i>sec~'  ( -^)  1  »  (5-8) 

obtained  after  substituting  (5-7)  into  (5-6) . 

From  (3-27)  the  maximum  azimuthal  mode  number  is  obtained 

v^,=V,  (5-9) 

in  the  limit  as  of-w. 
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This  section  is  provided  as  a  link  between  the  general 
formulation  of  Chapter  III  and  the  more  detailed  step  index 
analysis  contained  in  Appendix  E.  It  is  noted  that  (5-3)  ,  (5- 
7)  and  (5-8)  are  easily  shown  to  be  in  agreement  with  (E-20) , 
(E-21)  and  (E-23) ,  respectively. 

C.  ADDITIONAL  RESULTS 

The  results  and  discussion  appearing  in  Appendix  E  were 
presented  at  an  IEEE  conference  [Ref.  19]  .  Due  to  the 
limitation  of  space,  only  a  set  of  representative  trajectory 
curves  were  presented  at  the  conference.  The  parameters  for 
other  trajectory  examples  are  summarized  in  Table  7  in 
Appendix  E.  In  order  to  illustrate  the  application  of  the 
program  under  limiting  conditions,  two  additional  examples  are 
presented. 

1.  Example  A. 

As  seen  in  Figures  15  and  16,  the  rays  are  nearly 
meridional.  The  value  for  r^=0.0104,  which  for  m=30  and  t'=l, 

can  be  checked  by  use  of  (5-6)  .  It  is  noted  that  the  <f> 
dependence  on  r  is  given  by 

0=sec*'{^).  (5-10) 

■^min 

The  trajectory  lines  are  cibsolutely  straight  as  expected  in  a 
step  index  fiber.  In  Figure  15,  a  full  three-dimensional  view 
of  the  trajectory  is  provided.  It  is  noted  that  the 
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Figure  15:  Step  index  skew  rays  with 
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trajectories  continue  out  to  r=l,  which  can  be  seen  from 
Figure  16. 

2.  Example  B.  i'b191»i'„„b200 

Here  p  was  chosen  to  yield  a  value  r^  that  was  much 
closer  to  one,  i.e.,  a  high  rotational  component.  Again  for 
m=l  and  >'=191,  r„^=0.9358  can  be  checked  with  (5-6)  and  Figure 
17.  The  corresponding  three-dimensional  plot  is  Figure  18. 

D.  FINAL  REMARKS 

In  this  chapter  ray  trajectories  were  established  from  the 
general  analysis  of  Chapter  III  with  the  only  assumption  that 
From  Appendix  E,  we  can  see  that  the  model  relating  the 
mode  numbers  m  and  u  with  the  orientation  angles  6  and  J  is 
very  simple.  Also,  the  three-dimensional  perspective  of  ray 
trajectories  shows  that  the  model  provides  a  good  physical 
insight  for  a  deeper  understanding  of  optical  fibers. 
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VI.  NUMERICAL  SOLUTION  FOR  SKEW  RAYS  WITH  ARBITRARY  ot 

A.  INTRODUCTORY  REMARKS 

The  treatment  for  skew  rays  with  arbitrary  o;  requires  the 
full  generality  of  the  analysis  presented  in  Chapter  III. 
Meridional  rays,  which  are  covered  in  Chapter  IV,  will  not  be 
treated  here. 

B .  ANALYSIS 

For  the  power  profile  given  by  (3-13),  maximum  azimuthal 
and  radial  mode  numbers  can  be  predicted  from  (3-27)  and  (3- 

20)  ,  respectively.  Integral  limits  for  (3-20)  will  satisfy  (3- 

21)  . 

In  order  to  characterize  the  trajectories  given  by  (3-32) 
and  (3-34)  associated  with  the  p  and  m  modes,  the  normalized 
phase  constant  given  by  (3-17)  and  the  turning  points  need  to 
be  accurately  determined.  In  the  most  general  case  the  problem 
does  not  have  an  analytic  solution.  In  fact,  the  condition  on 
the  turning  points  r,  and  r^,  the  normalized  phase  constant 

e  (|8)  ,  and  the  radial  mode  number  given  by  (3-15)  need  to  be 
satisfied  in  a  self  consistent  manner.  The  conceptual 
flowchart  for  this  process  is  provided  in  Figure  19. 

Finally,  evolution  equations  can  be  predicted  from  the 
integral  form  of  (3-32)  and  (3-34) ,  which  become 


47 


Given  m^m 


max 


V  <V 
-  '^max 

®9e[0.1] 


1 

1 

1 

'  1 

Find  turning  points  (3-19) 

compute  m^ 

I 

i 

Compare  m  with  m 
6=  m-m^ 


Figure  19:  Calculation  of  e  (/3)  and  turning  points  for  skew 
rays  with  arbitrary  a. 
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-  _  yln^-NA^e  (0)  r? 

- 5a - j. 


rdr 


(6-1) 


e  (|8)r2-(^)2-(r)2- 


and 


> 


6(j8)r2-(^)2-r“*2 


dr 

r 


(6-2) 


respectively,  where  r€[r, .rj]. 


C.  EXAMPLES 

1.  Ray  trajectories 

Figure  20  shows  a  three-dimensional  plot  of  a  numerically 
calculated  skew  ray  trajectory  from  (6-1)  and  (6-2)  using  a=6. 
The  dense  vertical  lines  are  drawn  at  r=r, ,  and  the  low 

density  lines  are  drawn  at  r^r^.  From  Figure  20,  Az  is  defined 

as  the  axial  distance  between  ray  turning  points  r,  and  r2. 
Because  of  the  three-dimensional  perspective,  the  normalized 
increment  in  Az  only  appears  to  change  size.  Figure  21  shows 
the  corresponding  axial  view.  From  Figure  21,  A<p  is  defined  as 


A<^=2arctan  ( 


(6-3) 


Figure  22  shows  a  ray  trajectory  with  conditions  given  on 
the  Figure  with  As  can  be  seen  from  Figure  23,  which  is 
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Figure  20:  Ray  trajectory  for  conditions  as  stated. 
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the  axial  view  of  Figure  22,  for  this  case,  r,  is  close  to  rj 
as  expected  from  Figure  3  for  v  close  to 

For  the  ray  trajectories  shown  in  Figure  24  and  Figure  25, 
a  very  high  value  of  a  was  chosen.  As  can  be  seen,  the 
trajectories  are  indistinguishable  from  straight  lines,  as 
expected  for  rays  in  a  step  index  fiber  (as  discussed  in 
Chapter  V  for  step  index  case  in  the  limit  as  a-w) . 

The  equations  (6-1)  and  (6-2)  were  also  numerically 
integrated  for  the  case  a=2.  The  three-dimensional  perspective 
and  axial  views  are  shown  on  Figures  26  and  27,  respectively. 
It  is  noted  that  in  Figure  27,  it  is  clear  that  the  azimuthal 
orbit  is  closed.  These  results  will  be  reexamined  and 
explained  in  the  next  chapter. 

2.  Variation  of  fiber  parameters 

Table  2,  Table  3  and  Table  4  were  derived  from  repeated 
applications  of  the  optical  conditions  V=100,  NA=0.5,  n,=l.5 
and  q;  as  stated  on  the  tables.  The  corresponding  three- 
dimensional  representations  are  seen  in  Figures  26,  20  and  22, 
respectively.  The  trends  relating  and  Az  are  readily 
explained.  Because  all  of  these  tables  have  the  same  kind  of 
variation  of  parcimeters,  discussion  will  be  concentrated  on 
Table  2.  In  the  lower  half  of  the  table,  m  is  held  constant  at 
the  value  10  while  u  varies  between  4  and  16.  As  v  increases, 
the  angle  A<^  increases  with  the  rotational  component,  k^,  as 
expected.  It  is  clear  from  (2-13)  (or  in  normalized  form  from 
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Figure  26:  Ray  trajectory  tor  <x=2. 
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m 

V 

nW* 

^2 

Az 

A^trad] 

2 

20 

50 

2 

0.327 

0.610 

9.170 

0.984 

6 

20 

50 

4 

0.289 

0.689 

9.126 

0.794 

e 

20 

50 

6 

0.265 

0.754 

9.083 

0.675 

10 

20 

50 

8 

0.246 

0.812 

9.039 

0.588 

15 

20 

50 

15 

0.204 

0.978 

8.885 

0.411 

10 

4 

50 

23 

0.057 

0.690 

9.170 

0.164 

10 

8 

50 

21 

0.108 

0.740 

9.126 

0.289 

10 

12 

50 

19 

C  .152 

0.785 

9.083 

0.382 

10 

16 

50 

17 

0.193 

0.826 

9.039 

0.459 

Table  2;  Variations  of  fiber  parameters  for  ot-2 . 

(3-15)  and  (3-17)  that  with  m  fixed  any  increase  in  v  must  be 
compensated  by  a  decrease  in  |5,  The  interpretation  of  as  the 
z  component  of  the  propagation  vector  k(r)  (see  Figure  2)  then 
clarifies  the  slight  decrease  in  Az.  In  the  upper  half  of  the 
table,  V  is  held  fixed  while  m  varies  between  2  and  10. 


^max 

B 

nW 

m 

r, 

^2 

Az 

A(^  [rad] 

99 

2 

30 

10 

0.057 

0.999 

17.261 

0.114 

99 

10 

26 

10 

0.217 

0.999 

12.581 

0.437 

99 

18 

23 

10 

0.316 

0.999 

9.845 

0.643 

99 

20 

22 

5 

0.472 

0.999 

12.384 

0.984 

99 

20 

22 

15 

0.263 

0.999 

7.386 

0.532 

99 

20 

22 

20 

0.217 

1 

6.069 

0.437 

Table  3;  Variations  of  fiber  parameters  for  a=l0000. 
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V 

nW 

m 

r, 

^2 

Az 

A0 [rad] 

68 

5 

26 

5 

0.139 

0.705 

14.075 

0.396 

68 

30 

15 

5 

0.425 

0.850 

8.129 

1.030 

68 

45 

8 

5 

0.512 

0.902 

6.745 

1.176 

68 

5 

26 

25 

0.052 

0.986 

7.083 

0.105 

68 

5 

26 

15 

0.073 

0.879 

9.118 

0.166 

68 

5 

26 

5 

0.139 

0.705 

14.075 

0.396 

Table  4:  Variations  of  fiber  parameters  for  a=6. 

According  to  Figure  2,  it  is  clear  that  as  m  increases,  /S 
decreases,  and  as  such,  Az  should  decrease.  This  pattern  is 
again  evident  from  the  table.  Referring  back  to  the  relation 
between  /3  and  e  (|8)  ,  given  by  (3-17)  ,  and  the  discussion  of  the 
dependence  of  turning  points  on  e  (0)  in  Chapter  III,  it  is 
readily  confirmed  that  the  turning  points  will  move  together 
as  m  increases  (see  Figure  6) .  It  follows  from  (6-3)  that  if 
r,^r,+Ar,  and  rj-^rj-Ar^,  that  A<f>  will  decrease,  an  evident 

pattern  in  the  table. 

It  is  noted  that  Table  3  was  derived  under  the  saime 
conditions  as  the  step  index  skew  rays  for  Table  6  in  Appendix 
E,  except  that  here,  a  is  not  infinite  but  very  big.  The  fact 
that  the  results  of  those  two  tables  are  approximately  the 
same,  demonstrates  that  the  arbitrary  a  three  dimensional 
analysis  presented  in  this  chapter  is  in  agreement  with  the 
step  index  analysis. 
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D.  FINAL  REMARKS 


In  this  chapter,  ray  trajectories  were  established  for  the 
most  general  case  of  Chapter  III.  Although  many  numerical 
calculations  have  to  be  done,  the  model  is  still  fairly 
simple.  Also,  the  demonstration  that  trajectory  dependence 
the  mode  numbers  agrees  with  theoretical  expectations,  is  one 
more  sign  of  the  validity  of  the  model. 
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VII.  QUASI  ANALYTIC  SOLUTION  FOR  SKEW  RAYS  WITH  a»2 

A.  INTRODUCTORY  REMARKS 

In  this  chapter,  the  three-dimensional  skew  ray  analysis 
will  be  examined  for  the  special  case  a=2 .  This  case  has 
special  significance  for  telecommunications  applications.  It 
is  for  this  profile,  the  so  called  parabolic  profile,  that 
multimode  fibers  exhibit  an  approximate  minimum  in  modal 
dispersion  [Ref.  1].  It  turns  out  that  the  integrals 
describing  the  ray  trajectories  (see  (6-1)  and  (6-2))  and 
radial  mode  number  (see  (3-15)  with  f(r)=0)  can  be  evaluated 
analytically  for  the  case  a=*2.  Although  numerical  integration 
is  no  longer  needed,  the  general  relationship  between  m,  e  (fi) , 
r,  and  rj  is  too  complex  to  be  established  without  the 
assistance  of  a  numerical  root  finding  algorithm. 

B.  ANALYSIS 

For  Qf=2,  (3-27)  becomes 


The  radial  mode  number  can  be  found  by  setting  x=?  and 
f(r)=0  into  (3-15),  yielding 
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in= 


-T—l  €mx-{^)^-iedx. 


(7-2) 


where  the  turning  points  are  (see  (3-i9))  given  by 


€iff)± 


^1.2  =  - 


€^((8) -4(^)2 


(7-3) 


After  using  [Ref.  18] 


JVax^.bx^c  ^  =  - 


dx  dx  (7-4) 


!• 


\]ax^+hx+c  ^  xjax^+bx+c 


I 


dx 


1  2ax+b 


^ax^+bx+c 


;Sin'‘  ( . 


b^-Aac 


)  . 


(7-5) 


and 


/ 


^  -  =-l-sin-’(  ) 

x^ax^+bx+c  I  x|  ^b^-4ac 


(7-6) 


in  (7-2) ,  it  follows  that 


2Tr  2 


-2x+e  (|S) 


eH0)  -4(^)‘ 


.)  -.isin-'  (. 


e  {p)x-2[l) 


eM|8)-4(i:)2 


(7-7) 


where  the  limits  are  defined  by  (7-3) .  Equation  (7-7)  is  not 
analytically  invertible  so  it  is  necessary  to  have  the 
computer  find  e  (/3)  using  the  root  searching  algorithm  shown  in 
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Figure  19.  Similarly,  in^  can  be  found  after  setting  e(i8)=l  in 
(7-7)  and  directly  evaluating  the  expression. 

In  order  to  obtain  the  trajectory  for  a=2,  (6-1)  becomes 


'  ’  2jffl - k 


dx 


-x+e  {p)x-{l)^ 


(7-8) 


which  can  be  evaluated  using  (7-5)  and  leads  to 

_  \jn^-NA^e  (/S)  .  .  -2x+6(i8)  .,x 

- — - [-sin  ‘(__^=)]j.. 


2NA 


e"(i8)-(-^)' 


(7-9) 


where  x6[x,,X2]  .  The  limits  are  given  in  (7-3) 
For  <f){x)  ,  (6-2)  for  a=2  becomes 


dx 


X  -X 


x+e  (/3)x-(Z)^ 


(7-10) 


and  by  using  (7-6),  (7-10)  becomes 


€[0)X-2(L) 

<t>{x)  =-.i  [sin'*  (— 


eUP)  -4(-^)' 


:)]r. 


(7-11) 


where  x€[x, ,X2]  .  It  is  easy  to  show  that  substitution  of  (7-3) 
into  (7-11)  leads  to 


<#>(x,)  =0 


(7-12) 


and 
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(7-13) 


0(X2)=I, 

respectively.  After  four  applications  of  (7-11),  4>{x.)=2tt, 
which  proves  the  periodicity  of  ray  trajectories  with  qi=2,  in 
agreement  with  References  8  and  20. 

Figure  28  shows  a  three-dimensional  plot  of  a  numerically 
calculated  skew  ray  trajectory  for  conditions  given  on  the 
Figure.  Its  axial  view.  Figure  29,  demonstrates  clearly  that 
the  ray  trajectory  on  the  XY  plane  is  a  closed  orbit.  As  soon 
as  the  ray  trajectory  returns  to  the  same  point  on  the  XY 
plane  after  four  applications  of  (7-11),  the  z  component  of 
the  ray  trajectory  is  also  forced  to  be  periodic. 

Both  the  numerical  ray  tracing  program  used  in  Chapter  VI, 
and  the  quasi  analytic  ray  tracing  program  for  oi-2,  were 
tested  under  the  same  conditions  as  given  in  Figure  28.  It  was 
determined  that  both  programs  produced  exactly  the  same 
results  with  no  greater  then  0.1%  discrepancy.  These  results 
lend  confidence  to  the  reliability  of  the  models. 

C.  EXAMPLES 

Table  5  was  derived  from  repeated  applications  of  (7-9) 
and  (7-11)  under  optical  fiber  conditions  consistent  with  V, 
NA,  nj  and  a  defined  on  Figure  28.  The  parameter  Az  is  defined 
on  Figure  28,  and  A4>  is  defined  on  Figure  29. 

The  discussion  of  the  variation  of  parameters  in  Table  5 
has  already  been  explained  with  similar  tables  in  Chapter  VI. 
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Figxzre  29;  Ray  trajectory  for  oi=2. 
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m 

— 

nWx 

r, 

^2 

Az 

A(t)  [rad] 

2 

20 

50 

2 

0.327 

0.610 

9.170 

0.984 

6 

20 

50 

4 

0.289 

0.689 

9.126 

0.794 

8 

20 

50 

6 

0.265 

0.754 

9.083 

0.675 

10 

20 

50 

8 

0.246 

0.812 

9.039 

0.588 

15 

20 

50 

15 

0.204 

0.978 

8.885 

0.411 

10 

4 

50 

23 

0.057 

0.690 

9.170 

0.164 

10 

8 

50 

21 

0.108 

0.740 

9.126 

0.289 

10 

12 

50 

19 

0.152 

0.785 

9.083 

0.382 

10 

16 

50 

17 

0.193 

0.826 

9.039 

0.459 

5 

40 

50 

5 

0.447 

0.894 

8.885 

0.927 

Table  5:  Variations  of  fiber  parameters  for  a=2. 

It  is  easily  checked  that  the  parameter  variations  shown  in 
the  analytically  generated  Table  5  follow  the  same  trends  as 
seen  in  all  the  tables  of  Chapter  VI. 

D.  FINAL  REBSASKS 

In  the  previous  chapter  the  case  of  a=2  was  studied  with 
the  assistance  of  numerical  integrations.  It  was  demonstrated 
in  Figure  27  that  the  trajectory  actually  followed  a  closed 
orbit  for  the  specific  fiber  paraimeters.  It  was  shown  in  this 
chapter  that  the  result  is  exact  and  only  depends  on  the 
condition  that  a=2 .  This  somewhat  remarkable  result  has  been 
discussed  in  the  literature.  Specifically,  the  analogies 
between  mechanics  and  optics  provides  the  insight  as  the  basis 
for  this  effect  [Ref.  8] .  It  turns  out  that  only  two  kinds  of 
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potentials  exist  which  can  lead  to  closed  orbits.  First,  there 
is  the  ^  potential,  applicable  to  celestial  bodies,  and 

second,  the  r^  potential  associated  with  harmonic  oscillators 
[Ref.  20]  .  It  is  the  harmonic  oscillator  potential  which  is 
the  analog  of  the  parabolic  profile. 
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VIII.  CONCLUSIONS 


The  two  main  categories  that  deal  with  optical  fibers  are 
guided-wave  electromagnetic  modes  and  ray  optics.  Although 
these  approaches  both  describe  the  propagation  of  waves,  they 
are  not  directly  connected.  That  is  because  the  first  category 
relies  on  waveguide  theory  to  yield  mode  dispersion  fields; 
while  the  second  category  uses  initial  conditions  to  establish 
ray  trajectories.  Previously,  there  was  no  way  to  move  from 
angles  of  approach  to  waveguide  modes,  or  relate  ray 
trajectories  to  fields.  The  model  presented  in  Chapter  III  has 
established  a  general  formalism  for  connecting  waveguide  modes 
and  ray  trajectories.  Specific  cases  are  examined  in  Chapters 
IV  through  VII. 

As  shown  in  all  cases,  the  model  relating  the  mode  numbers 
is  fairly  simple.  This  and  the  fact  that  three-dimensional 
perspectives  of  ray  trajectories  can  be  easily  produced, 
facilitates  the  generation  of  physical  insight  for  a  deeper 
understanding  of  optical  fibers. 

It  is  noted  that  the  results  in  Chapter  V,  Chapter  VI  and 
Appendix  E  show  that  the  trajectory  dependence  on  the  mode 
nvimbers  agrees  with  theoretical  expectations.  The  direct  proof 
in  Appendix  B  that  shows  that  the  evolution  equations  are  in 
agreement  with  standard  Eikonal  analysis  makes  the  model 
presented  valid. 
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Finally,  the  connection  between  modes  and  rays  can 
ultimately  lead  into  new  methods  in  treating  waveguide 
dispersion  and  bending  losses,  but  this  is  a  matter  for 
further  research. 
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APPENDIX  A. 


CONCEPTUAL  FLOWCHART 


The  flowchart  shown  in  Figure  30  summarizes  the  main  steps 
in  the  development  leading  to  mode  specific  trajectories. 


Figure  30:  Conceptual  flowchart. 
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APPENDIX  B. 


A  COMPARISON  TO  EIKONAL  ANALYSIS 


The  proposed  evolution  equations  that  describe  the  ray 
trajectories  in  the  fiber  optic  according  to  the  Eikonal 
approach  [Ref.  9]  are 

(B-1) 

dz  E 

and 


^=_L 

dz  r^' 


(F-2) 


where  E  and  (  are  constants  of  the  motion  defined  as 


and 


B=n  ( r) 


dz 

di 


(B-3) 


(B-4) 

dz 

with  s  being  the  distance  of  the  point  measured  along  the  ray 
trajectory.  All  other  parameters  appearing  in  (B-1)  through 
(B-4) ,  are  defined  in  the  discussion  in  Chapter  III  and  Figure 
5.  In  the  discussion  to  follow,  we  demonstrate  the  relation 
between  the  Eikonal  constants  of  motion  and  the  invariants  0 
and  p . 

Definition  of  dl  in  Figure  5  and  (3-4)  leads  to 
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ds=dl . 


(B-5) 


By  substituting  (B-5)  and  (3-5)  into  (B-3) ,  we  have 


B=n(r) cos^ , 


(B-6) 


which  through  combinations  of  (2-10)  and  (3-2)  becomes 


(B-7) 


which  relates  the  phase  constant  3  to  the  Eikonal  constant  E. 
Also,  substitution  of  (3-5)  and  (3-7)  into  (B-4)  gives 


« _  sindsini 
^  ^  ■c5s"d“' 


(B-8) 


which  when  combined  with  (2-8),  (2-10),  (3-1)  and  (3-2)  leads 


(  =  -  " 

^  15- 


(B-9) 


Equation  (B-9)  relates  the  azimuthal  mode  number  v  and  the 
phase  constant  j8  to  the  Eikonal  constant  ( . 

Finally,  substitution  of  (B-7)  and  (B-9)  into  (B-1)  leads 


k^n 


(B-10) 


and  substitution  of  (B-9)  and  (B-10)  into  (B-2)  in  the  form 


dr/dz 

d4>tdz 


(B-ll) 


leads  to 
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p 

r 


_ 1 _ 

k^n^ir) 

r 


(B-12) 


Since  (B-10)  is  equivalent  to  (3-28)  ,  and  (B-12)  is 
equivalent  to  (3-33) ,  which  are  the  basic  evolution  equations 
of  the  proposed  analysis,  there  is  an  exact  agreement  between 
the  analysis  presented  and  the  Eikonal  analysis.  Nevertheless, 
the  treatment  presented  in  this  thesis  connects  Eikonal 
invariants  to  quantized  wave  mode  numbers. 
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APPENDIX  C. 


SYMBOL  TABLE 


NOTATION 

V 

a 

n, 

nj 

NA 

a 

k 


ki 


EXPLANATION 
V  parameter 
core  radius 
core  index 
cladding  index 
numerical  aperture 
power  profile  parameter 

vacuum  propagation  constant  { 

ton, 

core  propagation  constant  (  ) 


k2 


cladding  propagation  constant 


X 

p 

m 


r,0,  z 


e.i 


vacuum  propagation  wavelength  (2ir/k) 
azimuthal  mode  number 
radial  mode  number 

coordinate  system  defining  ray  location  (k, ) 
angles  which  define  ray  orientation  (k, ) 


e. 

/? 

6  (jS) 


critical  angle 
phase  constant 
normalized  phase  constant 
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r„un/  1^1.  ^2 

A 

Fioi) 

A 


extremes  in  radial  excursions 
trajectory  period 
normalization  integral 

index  difference  — 1 - 

2nf 


r,  z 


r  z 
a '  a 


i.  E 


Eikonal  invariants 
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APPENDIX  D.  PROPERTIES  OF  F (a) 


In  this  Appendix  we  examine  Fla)  given  by  (4-8),  that  is. 


Fla) 


=jVl-u“  du. 


(D-i) 


As  discussed  in  Chapter  IV,  Fla)  cannot  be  solved  analytically 
for  all  values  of  a.  However,  for  some  cases,  solutions  of 
F(a)  can  be  obtained  exactly. 

It  is  easy  to  show  that 


F(0)  =0 


(D-2) 


Floo)  =1, 


(D-3) 


Also,  by  using  [Ref.  18] 


fax+b 

3a 


(D-4) 


i^-x^  dx^iEVf^  +£!sin->  ( i5)  , 


(D-5) 


(D-1)  takes  the  form 


F{1)=| 


{D-6) 
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(D-7) 


F(2)  =1, 

4 

respectively. 

Since  then  for  a<l 

+  .  .=i.ainu, 

2  !  3  ! 

which  when  substituted  into  (D-1)  leads  to 

F,  ( or)  =JJv^aJnu  du . 

After  integration  of  {D-9)  [Ref.  18],  we  find 


(D-8) 


(D-9) 


(D-IO) 


for  Qf«l. 

For  a  much  greater  than  one,  and  since  |u(<l,  we  may 
assume  that 


(D-ll) 


Substitution  of  (D-ll)  into  (D-1)  and  integration  [Ref.  18] 
leads  to 


(D-12) 

2a+2 

for  c^l . 

Figure  8  shows  a  numerical  calculation  of  F(a)  for  several 
values  of  a.  Comparison  with  F,  (a)  and  F2(a)  shows  that  the 
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functions  are  asymptotic  to  F((x)  for  small  of  and  large  ot, 
respectively. 
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APPENDIX  E 


REFERENCE  19  AS  PRESENTED  AT  SSST  93 


A  Heuristic  Approach  to  the  Computation  of  3D-ray  Trajectories 

in  Step  Index  Optical  Fibers 


Athanasios  Nassopoulos  and  Ron  Pieper 


Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 
Monterey,  CA  93943 


Abstract 

3D-ray  trajectories  in  an  optical  fiber  are 
derived  through  use  of  a  simple  correspondence 
principle.  The  ray  trajectories  are  linked  to  the 
mode  number  of  the  exact  waveguide  solutions. 
The  analysis  and  simulations  presented  will  be  in 
terms  of  dimensionless  parameters  needed  to 
characterize  the  optical  fiber.  Specifically  this 
includes  the  V-parameter,  the  core  index  of 
refraction,  and  the  numerical  aperture.  The  curves 


prepared  are  presented  for  allowed  radial  and 
azimuthal  mode  numbers  and  are  presented  in  terms 
of  dimensions  which  are  normalized  by  the  core 
radius.  Although  this  approach  is  extendible  to 
graded  index  fibers  this  development  will  not  be 
presented  in  this  report. 

Introduction 

Exact  waveguide  solutions  have  been  derived  for  < 
select  number  of  index  of  refraction  profiles.  Amon( 
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these  exact  solutions,  the  derivation  for  the 
step-index  case  commonly  appears  in  most 
introductory  and  intermediate  level  optical  fiber 
texts  e.g.  [1].  Despite  the  mathematical 
importance  of  having  an  exact  solution  for  the 
propagation  in  an  optical  fiber,  the  wave  solutions 
can  be  difficult  to  interpret  physically.  It  is  not 
surprising  that  an  alternate  method  to  treating 
optical  fibers,  the  approximate  Eikonal  ray 
approach  [2],  which  in  turn  is  based  on  Fermat's 
extremum  principle,  has  been  repeatedly  applied  in 
order  to  visualize  effects  which  occur  in  optical 
fibers  [3].  For  example,  bending  losses  [41  and 
dispersion  effects  [5]  in  optical  fibers  have  been 
treated  using  this  approach.  However,  in  the 
Eikonal  approach  the  derived  ray  trajectories  are 
not  automatically  related  to  the  radial  and 
azimuthal  mode  numbers. 

Recently  an  approach  based  on  a  simple 
correspondence  between  waves  and  rays  has 
been  proposed  as  a  substitute  for  the  first  order 
WKB  analysis  in  producing  the  Eikonal-equivalent 
of  the  ray  propagation  equations  (6).  This 
correspondence  principle  has  been  used  to 
qualitatively  demonstrate  that  there  is  a  direct 


association  between  the  exact  wave  modes  and  the 
allowed  ray  trajectories  in  an  optical  fiber. 

Background 

Starting  from  the  Hemholtz  equation  which  is  the 
phasor  form  of  the  wave  equation  in  linear,  isotropic 
source-free  homogeneous  media,  we  have 

where  k(r)  is  the  wave  number  and  E  is  the  electric 
field.  It  is  noted  that  this  equation  is  only 
approximately  valid  in  a  graded  index  fiber. 
Consistent  with  the  form  given  in  (E-1 )  we  assumed 
that 

<^-2) 

where  k  is  the  vacuum  wave  number,  n(r)  is  the 
index  of  refraction,  and  r,  0  and  z  are  cylindrical 
coordinates  needed  to  describe  points  in  the  fiber. 
Following  a  standard  procedure  [1]  in  waveguide 
analysis  we  assume  that 

(E-3) 

where  v,  the  azimuthal  mode  number,  is  forced  to 
be  an  integer  due  to  periodic  boundary  conditions. 
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problem 

P  is  the  phase  constant  of  the  wave  guide.  F(r)  is 
the  radial  solution  to  (E-1|.  In  a  step  index  fiber 
the  exact  waveguide  solution  for  F(r)  can  be 
expressed  on  terms  of  Bessel  Functions  [1,3]. 

Using  a  simple  correspondence  rule,  V-^jk, 
it  is  easily  shown  [6]  that 


Equations  (E-4),  |E-5)  and  (E-6)  become  the  basi: 
for  the  three  dimensional  ray  modeling  that  will  b< 
the  focus  of  the  discussion  that  follows.  It  is  notec 
that  the  components  (k„  k,,  k,)  will  describe  thi 
direction  of  the  ray  and  that  (r,  0,  z)  locate  the  ray 
For  purposes  of  analysis  the  geometrica 
relationships  are  more  carefully  defined  in  the  nex 
section. 

Geometry  of  the  problem 

With  reference  to  Fig.  31  the  orientation  of  th< 
wave  vector  is  defined  in  terms  of  6  and  f.  I 
follows  that 


=  (E-4) 

*  r 

k,=p.  (E-5) 

A  condition  of  phase  synchronization  (11  on  the 
radial  part  of  the  wave  vector  results  in  an 
approximate  integral  expansion 


k^=k(i)9lx\Q  sin5  (E-7) 

k^=k(i)cosB  IE-8) 

k^=k(/)sine  cos^.  (E-9) 

Letting  dl  represent  the  line  variation  in  the  ra 
location  defined  by 


^rr, 


IE-6) 


dl=dz 


IE-10) 


where  the  radial  mode  number  m  takes  on  positive 
integer  values. 


and  that  the  incremental  components  d<p,  dr,  dz  ai 
tangemial  to  (E-7),  (E-8)  and  (E-9)  respectively,  it  follow 


03 


(/z»(//oos6  (E-11) 

dr=dl  sine  oosl[  2) 

r  sine  sln(.  IE-13) 

Tlirough  combinations  of  (E-7)-(E-13)  the  initial 
form  of  the  trajectory  equations  becomes 

dz  p  c/zie 
dr  k,  k,  oos5 


dr  k, 


r 


N 


=tanC. 


IE-15) 


The  evolution  equations  IE-14)  and  IE-15),  that 
describe  the  ray  trajectories,  can  be  shown  to  be 
equivalent  to  the  Ei)(onal  solution  [7].  Because  of 
the  3D  emphasis  in  this  report  we  have  chosen  to 
concentrate  on  skew  rays  1 1'  ^  0)  instead  of  the 
less  interesting  meridional  rays  I)*  -0). 


integration  IE-6)  must  remain  positive  would 
indicate  that  the  hatched  region  in  Fig.  32  defines 
the  allowed  extremities  of  the  possible  ray 
trajectories.  It  therefore  follows  that 


IE-16) 


Specification  of  a  radial  number  m  and  the 
azimuthal  mode  numbers  will  dictate  a  specific 
phase  constant  which  are  ordered  [1]  by 
convention 


where  the  upper  and  lower  bounds  on  are 
predictable  from  a  combination  of  wave  and  ray 
analysis  [11. 

The  turning  point  equation  for  r„in  is  obtained  by 
setting  the  integrand  in  IE-16)  to  zero 


IE-18) 

^Inln 


Modeling  Analysis.  Step  Index  case 

We  follow  the  usual  step  index  notation 
n|r)  =  ni  for  r^a  and  n|r)  =  n2  for  r>a.On  Fig.  32 
the  sum  of  the  two  terms  in  IE-6),  k*ni^-|»'/r)*,  is 
shown  plotted  versus  r.  The  condition  that  the 


The  maximum  allowed  v  is  determined  by  setting 
and  noting  from  Fig.  2  that  r„^„-*a.  It 
follows  that 

[ka^n^‘-n/]=lnt  (V^  IE-19) 

where  the  quantity  in  the  brackets  is  the  well 
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known  V  parameter.  The  tnt[  ]  operator  takes  the 
integer  part  to  produce  a  valid  mode  number.The 
general  expression  for  can  be  obtained  from 
(E-18)  as 


where  the  second  equation  follows  from  geometric 
considerations.  From  (E-5)  and  (E-8),  we  find  that 
(E-16)  can  be  expressed  as 

m(v)=//7f  [-  r  V(c  ) 

It  ■''ii*  r 

where  C  =  akniSin0,  r  =  r/a  and  r„i„  =  r„i„/a. 

integration  leads  to 

/77(v)=/nr  [l(v'c2-v2-vsec-’(-))]. 
n  V 


To  find  m^{i»)  it  is  necessary  to  allow  C  to  take  its 
maximum  value.  As  seen  from  (E-5)  and  (E-8)  this 
occurs  for  the  minimum  It  follows  from  (E-17) 
that  the  maximum  value  of  C  is  the  V  parameter. 
After  direct  substitution  into  (E-22),  m^(i^) 
satisfies 

It  V 

where  for  p-0 

/7w(0)=/n<-^l  <E-24) 

jr 

in  agreement  with  a  mode  space  analysis  [8]. 

From  the  set  of  equations  (E-20)  and  (E-22)  it  can 
easily  be  shown  that 

nTii=v[  (t— )*-1 -sec'’(T— )1  (E-251 

Given  that  m  and  v  are  known,  the  exact  value  of 
can  be  found  by  numerically  solving  (E-25). 

Normalized  ray  trajectories 

In  this  section  the  normalized  trajectories  for  rays 
defined  by  a  combination  of  standard  optical 
parameters  (V,  NA,  n^)  and  mode  indices  m  and  v. 
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Figure  33:  3D  representation  showing  4  turning  points 


Consistent  with  the  purpose  of  this  report  the 
results  for  meridional  rays  ( >>  =  0)  will  be  presented 
without  significant  discussion. 

Starting  from  IE-41  and  IE-7)  it  easy  to 
demonstrate  that 

y>  =  -r  k^=-r  k  n.,  sine  sin^.  IE-26) 

After  substitution  of  IE-20)  into  IE-26)  we  have 
that 

sin5=-^,  <E-27) 

r 


The  integral  form  of  IE-15)  which  is 

=  tanErfr= / "  — ^  (E-28) 

can  be  solved  by  substituting  IE-27).  Finally  we 
have 

<t.W=sec-’(3^)  IE-29) 

^inln 

which  is  not  valid  for  v=0.  However  from  IE-15) 
and  IE-27)  0|r)  =  O  for  j'=0.  z|r)  can  be  found  by 
solving  the  integral  form  oflE-14)  which  is 


which  is  needed  in  the  analysis  to  follow. 


Figure  34:  Projection  in  polar 
plane . 


^h: 


ctnd 

oos( 


dr. 


(E-30) 


With  a  second  application  of(E-27)  we  find 


where  z  «z/a.  The  angle  0  can  be  easily  calculated 
from  IE-20)  as 


sine= 


(E-32) 


and  where  r„i„  can  be  obtained  from  (E-25).  For 
i'=0  it  can  be  shown  that  z(r)=fctnS  and 
sin0  =  (mNA)/|Vni)  are  applicable  to  trajectory 
calculations  of  meridional  rays. 

Along  with  the  input  parameters  V,  NA,  n,, 
p  (given  from  (19)  )  and  m  (given  m^  from 
(E-23) ),  the  set  of  equations  (E-25),  (E-29),  (E-31 ) 
and  (E-32)  determine  the  normalized  ray 
trajectories  of  a  skew  ray  in  a  step  index  fiber. 


Discussion  of  results 

Figure  33  shows  a  3D  plot  of  a  numerically 
calculated  skew  ray  trajectory  for  conditions  given 
on  the  figure.  The  dense  vertical  lines  are  drawn  at 
r=7„„„  and  the  low  density  lines  are  drawn  at 
7=  1 .0.  Points  1 ,  2,  3,  4  are  turning  points  at 
7=1.0.  Because  of  the  3D  perspective  the 
normalized  increment  in  Az  appears  to  change  size. 
Figure  34  shows  the  corresponding  axial  view.  The 
angle  is  defined  by  rays  joined  at  the  turning 
points. 

Table  6  was  derived  from  repeated  applications 
under  optical  fiber  conditions  consistent  with  V,  NA 
and  n,  defined  on  Figure  31.  The  trends  relating 
variation  of  At^  and  Az  with  p  and  m  are  readily 
explained.  In  the  upper  half  of  the  table,  for  fixed  m, 
the  angle  A(p  increases  with  the  rotational 
component,  k^,  as  expected.  Also,  since  m  (E-6)  is 
constant,  0  must  decrease.  This  is  manifest  by  Az 
decreasing  with  p.  In  the  lower  half  of  the  table, 
where  p  is  fixed,  as  m  increases,  0,  as  seen  from 
(E-6)  must  again  drop.  Again  Az  must  therefore 
decrease.  And  finally,  the  parallel  decrease  in  L0  (E- 
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At^fradl 

Az 

m 

9 

0.12 

17.25 

10 

02 

0.44 

12.57 

10 

10 

0.64 

09.84 

10 

18 

0.98 

12.37 

05 

20 

0.53 

07.38 

15 

20 

0.44 

06.10 

20 

20 

Table  6  .-Representative  variation 
in  ray  parameters  with  wave  mode 
numbers . 

29)  can  be  explained  geometrically  as  seen  from 
figure  34  by  a  required  decrease  in  r^j^.  This  is 
apparent  by  noting  that  as  the  radial  mode  number 
increases  with  v  fixed  the  phase  constant  p 
should,  as  seen  from  (E-61  decrease.  Conditions 
(E-5)  and  (E-8)  then  indicate  6  will  increase  and 
therefore  from  (E-20)  will  decrease. 

Conclusions 

The  main  point  of  this  report  on  optical 
fiber  ray  calculations  is  to  demonstrate  the 
existence  of  a  pedagogically  attractive  alternative 
to  more  formal  methods  such  as  WKB  or  Eikonal 


analysis.  An  important  feature  of  the  method 
presented  is  that  it  links  the  traditional  cylindrical 
waveguide  modenumbers  for  the  exact  solution  to 
the  ray  trajectories.  Because  of  this  link  the 
numerical  implementation  of  the  model  can  provide 
valuable  physical  insight  for  the  fiber  optic 
application  areas  such  as  sensors  or  lightwave 
telecommunications. 
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